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ABSTRACT:  Coronal  images  contain  many  complex,  superposed,  faint, 

optically-thin  structures:  rays,  loops,  curved  streamers,  etc.  The 
spatial  resolution  is  often  limited  by  the  signal/noise  ratio  in  the 
recording  medium  while  the  dynamic  range  in  intensity  is  severly 
dominated  by  the  radial  gradient.  We  present  the  results  obtained  with 
a  new  simple  image  processing  algorithm  applied  to  several  typical 
eclipse  and  coronographic  pictures.  This  2D  filter  (second  order  diff¬ 
erence  operator  in  k-di section)  shows  superior  quality  by  improving  the 
signal/noise  ratio,  reducing  the  dynamic  range  and  enhancing  faint 
coronal  structures.  A  comparison  with  results  obtained  with  more 
commonly  used  2D  filters  is  also  presented.  We  describe  briefly  the 
algorithm  used  here,  together  with  discussions  of  enhanced  coronal 
images  of  coronal  physics  interest. 
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1.  Introduction 


The  extraction  of  precise  parameters  (intensity  distribution  over  a 
cross-section,  shape,  size,  direction,  etc.)  from  coronal  images,  both 
for  white-light  and  emission-1 ine  studies,  is  difficult  unless  some 
image  processing  technique  is  applied.  For  example, the  structures  to  be 
studied  are  optically  thin  and  superposed  structures  can  reproduce 

rather  large  modulations.  Additionally,  the  steep  radial  gradient  of 
intensities  can  mask  details  of  the  inner  coronal  structure,  while  other 
details  can  be  lost.  Further,  because  the  radiations  are  faint  and 

because  the  most  popular  detector  in  coronal  imaging  is  photographic 

film  with  its  rather  low  quantum  efficiency,  (the  use  of  CCD  cameras, 
with  all  their  special  characteristics,  is  not  considered  here), 
pictures  are  often  noisy,  with  variable  si gnal -to-noise  ratio.  Here  we 
present  the  results  of  some  attempts,  including  the  use  of  an  apparently 
novel  algorithm,  to  improve  coronal  pictures  in  order  to  try  to  answer 
several  outstanding  questions  concerned  with  different  aspects  of  the 
coronal  physics. 

Even  in  the  case  of  images  recorded  in  the  brightest  visible  line, 
530. 3nm  of  Fe  XIV,  some  image  processing  can  be  valuable;  the 

interpretation  of  detailed  structure,  such  as  interacting  post- flare 
loops,  is  otherwise  extremely  difficult,  and  it  is  virtually  essential 
in  comparative  line  studies,  especially  in  investigations  of  temperature 
and  density  inhomogeneities.  Figure  1  shows  a  part  of  a  picture  taken 
by  one  of  us  (J.S.)  during  the  solar  total  eclipse  of  July  31,  1981, 
with  a  3A  bandwidth  530.3  nm  filter.  The  picture  was  processed  using 
different  methods,  including  the  commonly  used  method  of  high  pass 
filtering  (see  figure  2).  Based  on  the  results  of  applying  such 

different  methods  of  processing,  largely  described  in  books  like  the 
classical  monograph  of  Rosenfeld  and  Kak,  1976,  or  those  of  Coster  and 
Chermant,  1985,  we  have  investigated  a  more  empirical  approach,  taking 
into  account  the  anisotropic  properties  (due  to  the  dominance  of 
magnetic  field  effects  and  the  dynamical  characteristics  of  the  corona) 
of  the  structures  in  the  solar  corona. 

An  example  of  the  result  obtained  with  this  apparently  novel  algorithm, 
which  we  call  "MaD  MAX  II"  (multi -directional  maximum  of  second 

derivatives)  is  shown  on  figure  1,  part  d,  as  compared  with  the  results 
of  applications  of  more  commonly  used  algorithms.  Since  we  believe  our 

algorithm  is  superior,  we  present  a  brief  description  of  the  method 

before  commenting  on  its  application  in  several  interesting  cases,  as 
well  as  some  unsuccessful  attempts. 


ure  1.  Example  of  image  processing  of  a  tvpical  coronal 
picture  using  different  algorithms: 

a.  Original  picture  taken  during  the  July  31st,  1981 
Solar  eclipse  (  J.  Sykora  )  in  the  light  of  Fe  XIV 
emissions  at  530.3  nm; 

b.  Image  obtained  with  a  laplacian  filter; 

c.  Image  obtained  using  the  "multidirectional 
maximum  of  absolute  derivatives" 

(J.M.  Malherbe  and  al.,1986  ); 

d.  Image  obtained  using  the  proposed  algorithm 
(  madtnax  ). 

ihat  mgu  pass  or  p»ss  bund  fillers  give  more 

FOCIlltc*  O 
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2-  Brief  description  of  the  proposed  image  processing 


The  purpose  of  the  processing  is  to  find  ,as  far  as  possible,  the  detailed 
distribution  of  different  features  (curved  or  straight,aligned  or  not)  of 
variable  sizes  which  typically  intersect.  These  features  axe  distributed 
over  a  noisy  picture  with  a  variable  in  intensity  background. For  such  an 
application, the  use  of  the  second  derivatives  technique  seems  appropri¬ 
ate.  For  digital  pictures  differences  must  be  used  rather  than  deriva¬ 
tives.  The  features  can  be  oriented  in  any  direction;to  find  this  orienta¬ 
tion  we  compute  the  second  order  differences  in  8  directions  and  choose 
those  for  which  the  values  are  maximum. 

We  shall  now  show  in  more  detail  how  the  processing  works. 

Let  g(x,y)  be  the  matrix  of  measured  intensities.  The  algorithm  uses 
the  following  operator  M  defined  by: 

M(g)=aM8(g)+^E2(M8(S2(g)))+'yE4(M8(S2(S2(g))))  +  ... 

where: 

->  M8  =  max(  -  A*1  ) 

with  A  **  the  second  order  difference  operator 

in  the  k  -direction  defined  by  the  angle  k  II/8,  k=0,l,...,7. 

->  S2  is  the  operator  which  averages  every  2  row*  and  every  2  columns 
of  a  matrix  to  make  a  single  row  and  a  single  column. 

->  E2  and  E4  are  the  operators  which  expand  a  matrix  by  bilinear 
interpolation  with  the  row  and  column  expansion  factor  equal  to  2  and  4 
respectively. 

->  a  ,  /?  ,  “7  are  some  weights  taking  the  values  0,1,2,...  depending  on 
the  signal  to  noise  ratio  and  on  the  relative  dimensions  of  the  picture 
features. 

The  operator  M  is  non  linear,  it  is  based  upon  the  anisotropic  opera¬ 
tor  max  and  it  is  shift-invariant  and  non  idempotant. 

The  effect  of  M  on  a  picture  can  be  understood  by  considering  a  sim¬ 
ple  one-dimensional  example.  Then  M8  is  computed  only  in  one  direc¬ 
tion  and  its  effect  will  be  to  amplify  "convex"  and  "concave"  points. 
Thus  its  purpose  is  the  detection  of  places  where  there  is  a  significant 
change  in  gray  level.  So  the  meaning  of  the  processing  in  terms  of  image 
sharpening  is  quite  obvious. 

Unfortunately, the  algorithm  uses  a  second  order  difference  operator 
relatively  sensitive  to  image  noise,  although  the  averaging  operator  S2 
reduces  this  effect. 
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Figure  2.  Image  processing  with  a  commonly  used  2D  filter 
(high  pass)  of  2  simultaneous  eclipse  monochromatic 
pictures  obtained  on  July  31st,  1981  (  J.  Sykora  ): 

a.  Original  image  in  637.4  nm  ofFe  X 
(  FWHM  of  the  filter  0.3  nm): 

North  is  At  top  left.  East  is  at  bottom  leftj 

b.  Processed  image  of  a.; 

c.  Original  image  in  530.3  nm  of  Fe  XIV 
(FWHM  of  the  filter  0.2  nm  ); 

d. rrocessed  image  of  c.; 

Note  the  absence  of  correlations  between 
the  emissions  of  Fe  X  (  T  10  K  )  and 
the  emissions  of  Fe  XIV  (  T  ~1.8xl0  K  ). 


3.  Discussion  of  results 


a.  Figure  2  shows  full  images  obtained  simul taneously  in  the  green  line 
(530. 3nm)  Fe  XIV  emitting  preferentially  at  temperatures  near  2.0  x 
106  K  (part  c)  and  in  the  red  line  (637. 4nm)  of  Fe  X  emitting  near 
1.2  x  10°  K.  Processed  images  (part  b  and  d)  with  a  "classical" 
high-pass  filter  in  the  Fourier  space  removes  the  radial  gradient 
and  shows  convincingly  the  large  spatial  differentiation  of 
locations  emitting  these  two  lines. 

This  makes  rather  doubtful  results  based  on  the  use  of  the  so- 
called  line-ratio  method  to  deduce  the  coronal  temperature  and 
abundances:  structures  emitting  in  one  temperature  "regime"  do  not 
spatially  coincide  with  structures  emitting  in  another,  when 
ionization  temperatures  differ  by  a  factor  of  1.5  or  more,  as  is  the 
case  for  Fe  X  and  Fe  XIV. 

Using  the  MaD  MAX  II  program,  we  succeeded  in  showing  the 
correlations  between  530. 3nm  emission  of  Fe  XIV  and  structures 
revealed  by  a  corresponding  white  light  picture  (essentially 
electron  density  Ne,  structures , ) .  Figure  3  shows  a  part  of  the 
530. 3nm  picture  processed  with  MaD  MAX  II,  as  was  the  W.L. 
picture.  Although  the  intensity  modulations  of  the  original 
pictures  are  very  different  (530. 3nm  emissions  are  a  nearly 
quadratic  function  of  the  Ne  distribution,  while  W.L.  intensities 
are  only  a  linear  function  of  Ne).  In  addition,  the  two  pictures 
were  obtained  using  completely  different  methods  and  instruments. 
The  structures  revealed  by  the  processing  look  very  similar, 
confirming  the  excellence  of  the  Fe  XIV  line  as  an  indicator  of  the 
coronal  density  distribution  and  therefore,  as  a  "good"  line  to  be 
used  for  monitoring  coronal  holes  (see  R.  Altrock,  this  volume).  We 
do  not  know  the  exact  origin  of  this  correlation,  but  suspect  that 
the  most  evident  interpretation  calls  for  a  temperature  "regime"  of 
the  solar  corona  that  is  optimum  for  Fe  XIV  emission,  although  this 
interpretatin  should  be  tested  using  several  other  coronal  lines 
such  as  those  of  Fe  XIII  in  the  UV  or  IR. 

b.  The  use  of  the  proposed  algorithm  appears  to  be  very  effective  in 
the  case  of  W.L.  pictures,  showing  subtle  details  connected  with 
dynamical  processes  such  as  coronal  mass  ejections,  etc.  Figures  4 
and  5  show  the  images  processed  with  MaD  MAX  II  of  a  coronal  mass 
ejection  phenomena  observed  during  the  Feb.  16,  1980  total  solar 
eclipse  (the  so-called  "tenni s-racket"event) ;  this  CME  is  studied 
elsewhere  in  this  volume  (see  Koutchmy  and  Nitschelm),  so  we  do  not 
comment  here  on  the  scientific  aspects  of  the  result  of  the  image 
processing.  Observe,  however,  that  the  unprocessed  picture  does  not 
reveal  the  event  convincingly  (compare  with  figure  5).  However,  the 
processed  images  (Figures  4  and  5)  suggest  an  interpretation  based 
on  the  twisting  or  reconnection  processes  and  facilitate 
measurements  of  proper  motions. 

c.  The  last  example,  Figure  6,  is  an  example  of  very  interesting  photo¬ 
graphic  images  of  Fe  XIV,  FeX  and  Ha  emissions  obtained  with  the 
"one-shot"  coronagraph  of  Sacramento  Peak  (NSO)  Observatory  (see 
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Figure  3.  Comparison  of  different  coronal  emissions  of 
the  inner  corona: 

Upper  part:  processed  image  of  530.3  emissions  of 
Fe  XIV; 

Lower  part:  processed  image  of  white  light  intensities 
(density  structures)  using  different  weights  (a&b). 
Note  the  excellent  correlations  between  the  re  XIV 
emissions  and  the  density  structures  as  shown  on 
the  white  light  picture. 
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I  igure  4.  Nlorphologicul  analysis  of  the  legs  of  a  coronal  inass 
ejection  event  observed  during  the  total  solar  eclipse 
of  f  elt.  1  0. 1  9 HO  (  see  Fig, 4  ) : 

a.  Part  of  an  original  picture  taken  by  J.  Diirst 
at  10b 1 0  IT. 

lj..c..d.  1'rocessed  images  with  different  weights  in 
the  proposed  algorithm  (  "madmnx"  ). 

Note  the  twistings  apparent  in  the  processed  images  at 
the  location  of  the  legs  of  the  O.M.K.  shown  on 
t  he  Fin.'*. 


Figure  5.  Morphological  analysis  of  C.M.E.  observed  during 
the  total  solar  eclipse  of  Feb. 10,  1980: 

a.  Part  of  an  original  picture  taken  by  F.  Diego 
at  8 h 2 .'5  U.T.  _ 

b. ,c.  Processed  images  with  different  weights  in 
the  proposed  algorithm  (  "madmax"  ). 
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Smartt  and  Zhang,  1986),  which  have  not  yet  been  successfully  proc¬ 
essed.  These  pictures  are  recorded  on  2415  film  of  very  fine  grain 
with  a  scale  of  55  arcsec/mm.  For  several  reasons,  it  nas  not  been 
possible  to  digitize  these  images  successfully  with  the  fast 
microphotometer  of  SPO.  Primarily,  the  small  scale  of  the  images  is 
not  fully  resolved.  An  additional  difficulty  comes  from  the  so- 
called  "fringing  effect"  produced  by  the  laser  light  used  to  scan 
the  picture.  Further,  the  overall  density  of  the  film  is  low  (d< 
0.2)  and  the  range  of  photographic  densities  to  be  measured  are 
small.  We  applied  MaDMAXII  to  several  digitized  pictures  and  even 
more  "classical"  algorithms  without  great  success.  Disappointing 
results  were  obtained  even  when  a  magnified  duplicate  was  used  for 
the  digitization.  Figure  6  shows  selected,  highly  magnified  prints 
of  images  of  great  interest  (other  very  interesting  examples  have 
also  been  studied,  but  they  represent  only  a  very  small  fraction  of 
the  data  set  available  from  this  instrument.  The  sketch  illustrates 
the  kind  of  events  that  can  only  be  studied  accurately  with 
processed  images.  It  shows  Ha  emission  of  the  cool  part  of  the 
corona  (indeed  -  1 oop-1 ike  prominence  structures)  contai ned  within 
the  system  of  coronal  loops  as  evidence  by  their  emissions  in  Fe  XIV 
line.  Such  observations  strongly  suggest  a  "2-step"  configuration 
of  the  magnetic  field  lines  above  an  active  region,  substantially 
supporting  the  phenomenology  and  mechanism  of  CMEs  recently 
suggested  to  explain  the  observation  of  prominence  material  behind 
transients.  However,  this  somewhat  rather  speculative  claim  needs 
to  be  checked  by  comparing  the  morphological  character!' sti cs  of  the 
structure  as  revealed  in  the  processed  images.  But  for  this  we  need 
an  improved  method  for  obtaining  accurate  digitized  data. 
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